Time evolution of the QED vacuum in a uniform electric field: 
Complete analytic solution by spinorial decomposition 
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Exact analytical solutions are presented for the time evolution of the density of pairs produced in 
the QED vacuum by a time-independent, uniform electric field. The mathematical tool used here 
to describe the pair production is the Dirac-Heisenberg-Wigner function introduced before [Phys. 
Rev. D 44, 1825 (1991)]. The initial value problem for this function is solved by decomposing the 
solution into a product of spinors. The equations for spinors are much simpler and are solved analyt- 
ically. These calculations are nonperturbative since pair production is due to quantum-mechanical 
tunneling and the explicit solutions clearly exhibit their nonanalytic behavior. 



I. INTRODUCTION 



The problem of pair creation from the QED vacuum by 
electromagnetic fields has been studied by various meth- 
ods since the early days of quantum electrodynamics. 
The application of the one-particle Dirac equation leads 
to difhculties (Klein paradox) because pair creation is ob- 
viously a many-body problem. All paradoxes are resolved 
when this process is described by quantum field theory. 
Since pair creation by an external electromagnetic field 
involves quantum-mechanical tunneling, the process can- 
not be described in the framework of perturbation theory. 
Some time ago we introduced in [l| a formalism based 
on a generalization of the Wigner function to quantum 
field theory. This allowed us to describe pair creation 
without resorting to perturbation theory. The basic tool 
in this formalism is the 16-component Dirac-Heisenberg- 
Wigner (DHW) function. The formalism based on the 
DHW function was further developed in j^-Q to yield 
quantum transport theor y fo r relativistic spinning parti- 
cles. In several papers [6l-[ll| a relativistic generalization 
of the Wigner function was introduced. However, this 
approach is not suited for the solution of the initial value 
problem which is the main subject of this work. Other 
methods were also developed to treat the problem of pair 
production p^ - [T^ but they have not led to an exact an- 
alytic solution. Such a solution is especially useful since 
the dependence on the field strength of the pair density 
is described by a function with an essential singularity. 

In the general case of space-dependent fields the evolu- 
tion equations for the components of the DHW function 
form a set of complicated integro-differential equations. 
However, for a time-independent, uniform electric field 
these equations become tractable since, as shown in [l[, 
they can be reduced to the following three equations for 
the relevant components of the DHW function (notation 



is explained in the next section): 
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These equations were integrated numerically in [l| and 
the Schwinger formula [Tsj for pair creation was obtained. 
Recently, by an elaborate procedure a single solution of 
equations related to Eqs. ([ij was constructed as a combi- 
nation of parabolic cylinder functions [l^ . This discovery 
stimulated our return to the study of pair production by 
a uniform electric field. 

In the present paper we construct not one but a com- 
plete set of solutions of Eqs. (P). With the use of these 
solutions we solve the initial value problem for these 
equations with arbitrary initial conditions. The initial 
conditions that interest us most are those of the QED 
vacuum and we describe the time evolution in this case. 
Our solution clarifies the mechanism responsible for pair 
production and may explain some details of the vacuum 
decay studied recently |17| . 

In Sec. |TT] we recall basic definitions and we derive a 
simplified form of the evolution equations (HJ best suited 
for our analysis. In Sec. IIIII we exhibit the underly- 
ing 0(3) symmetry that enabled us to express the solu- 
tions of the evolution equations in terms of spinors. The 
method of spinorial decomposition introduced in Sec. IIVI 
to lower the order of the differential equations by reduc- 
ing the number of functions is quite general and might 
be applicable to a class of similar problems. Finally, in 
Sec.|V]we describe the time evolution of the QED vacuum 
and we give physical interpretation of our results. 



II. THE DHW FUNCTION AND ITS 
EVOLUTION 
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field was defined in [l| as follows: 



10 matrix M (p) has the following block structure 



*a {r + s/2,t),^ ir~s/2,t) 



(2) 



where the square bracket denotes the commutator of the 
quantized Dirac field operator 'S' with its hermitian con- 
jugate and the phase if is the line integral of the vector 
potential 
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where we introduced the following shorthand notation for 
the 3x3 matrices: 



The commutator was introduced to make the DHW func- 
tion invariant under charge conjugation. The phase fac- 
tor ensures its gauge invariance and, as a result, the mo- 
mentum variable p is the kinetic momentum and not the 
gauge-dependent canonical momentum. The DHW ma- 
trix Wa/3 will be expanded into a set of Dirac matrices, 
as was done in [l| 
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W{r,p,t) - - I l4/o(r,p,t)+^p,/,(r,p,t) 



-(- <T-go (r,p,t) -f- ^Pjcr-gj (r,p,t) ] . (4) 

i=i 



Since the DHW matrix is hermitian, the expansion coef- 
ficients are all real functions. The Dirac equations for the 
field operators lead to evolution equations for the DHW 
function W that can be rewritten as 16 equations for the 
functions fi and gi (i = 0, . . . , 3). 

For arbitrary electric E{r,t) and magnetic B{r,t) 
fields the evolution equations were derived in 1] . Here we 
shall need only a simple version of these equations valid 
for a time-independent, uniform electric field. Since we 
will be concerned only with the time evolution of the ini- 
tial vacuum state, it is sufficient to consider only ten, 
out of 16, components of the DHW function that cou- 
ple to the vacuum DHW function. The relevant com- 
ponents can be combined into a ten-dimensional vector 
W = (/3 , 9o , 9i , 92 ) and the evolution equations have the 
form: 



{dt + eE-dp) W {p, t) = jM (p) W (p, t) , (5) 



where we have taken into account that a uniform electric 
field does not introduce any space dependence. The 10 x 



/,From now on we shall use units c = I = h. 

The vacuum DHW function contains only the compo- 
nents = —2m/Ep and gi = —2p/Ep, where Ep = 
^ym^ +p^. The solution of the evolution equations that 
begins as the vacuum state, as shown in [l[, lies in a 
three-dimensional subspace and it can be described by 
the three-dimensional vector (wi(p, t), ■W2{p, t),W3{p, t)), 

W {p,t) = -2 [wi{p,t)ui + W2{p,t)u2 + W3{p,t)u3] , 

(7) 

where the ten-dimensional orthonormal vectors Uk , in the 
same block convention as used for M{p), have the form 

ui = — [ m I I Pi I ] , (8a) 

m± 
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[ I p± X n I I mn ] , 



ii3 = [0|0|n|0]. 



(8b) 
(8c) 



where m± = \J rr? -I- ■p\ is the transverse mass (energy) , 
n is the unit vector in the field direction, and bold ze- 
ros stand for three-dimensional zero vectors. Note that 
the derivative in the field direction present in ([5]) does 
not affect the vectors since they depend only on the 
transverse momentum p^. The DHW function of the 
vacuum state in this basis is 



2to^ 2p|| 

Wvac = "1 - -ET^s, 

£jp £jp 



(9) 



where : 



np. The subspace spanned by the vectors Uk 
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is invariant under the action of the matrix M(p), namely 

M(p)wi -p||M2, (10a) 
M{p)u2 = p\\ui - TOj^ua, (10b) 
M(p)m3 = m_LU2- (10c) 

With the use of these relations, by comparing the coeffi- 
cients in Eq. ([5]) that multiply the vectors Uk^ we obtain 
the following three equations for the coefficient functions 



(11) 



where we introduced the following dimensionless vari- 
ables (c and h are reintroduced below for clarity) 
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These evolution equations are equivalent to the original 
ones JT]), treated in [l[ but their new form is more conve- 
nient. The new functions (wi,W2,u'3) are related to the 
old ones (00,01,02) by a rotation. 



Op - goi 



W2 = 02, W3 



gap + ai 



(13) 



The solution of the evolution equations pT|) will be re- 
duced in the next section to the solution of ordinary dif- 
ferential equations. 



III. REDUCTION TO A PRECESSION 
EQUATION 

The solution of the initial value problem for the equa- 
tions (ITT|) will be reduced now to a set of ordinary differ- 
ential equations by an extension of the method of char- 
acteristics. To this end let us denote by Xk{q) three 
linearly independent vectors that are solutions of the fol- 
lowing set of equations 



f-Arfc(g) = L{q)Xk{q), 



(14) 



where L{q) is the matrix appearing in Eq. (jlip and k 
labeling different solutions runs from 1 to 3. From the 
vectors Xk{q) we can build a nonsingular matrix W(q), 



W('Z) 



vi vi vi 

^1 '^2 '^3, 

v2 v2 v2 

'^l ^2 ^3 

>'3 v3 v3 

^1 ^2 ^3 



(15) 



where the lower index labels different vectors and the 
upper index labels the vector components. The solution 
of the initial value problem for Eq. pT|) has the form 

X{q, t) = W(g) \w{q - Et)] X{q - Er, 0). (16) 



The validity of this formula can be checked by inspec- 
tion. First, we check that this expression is a solu- 
tion of Eq. ()11|) . The combination of partial derivatives 
{d-r -f Edq) when acting on the first part of reduces 
to the derivative with respect to q while it does not af- 
fect the part dependent on q — Et. The derivative with 
respect to q reproduces the right hand side of equations. 
Second, we check that this expression satisfies the initial 
condition. This follows from the fact that at t = both 
arguments become equal and the product of the matrix 
and its inverse reduces to the unit matrix. 

To make use of the formula (jl6p we must find a com- 
plete set of solutions of Eq. (IT4l) . The reduction of the set 
of three equations ([T4| to just one equation for a single 
function leads to a third-order equation. In contrast to 
second-order equations, such equations are, in general, 
hard to solve. A representation of the solution of the 
resulting third-order equation in the form of a contour 
integral has been given by Radozycki (IS] . The contour 
integral form was sufficient to obtain the Schwinger for- 
mula. However, it cannot serve as a basis for the solution 
of the initial value problem and we proceed by a different 
route. 

The key property is the presence of the 0(3) symme- 
try in the equation Indeed, this is an equation de- 
scribing the precession in three-dimensional space of the 
vector X{q) around the vector 9i{q), 



EX'{q) = -2n{q)xX{q), 



(17) 



where H = {1,0, q) and q plays the role of time. The 
factor of 2 is pulled out for further convenience. Were it 
not for the g-dependence of the precession vector 'H{q), 
the solution of Eq. ((T7)) would have been very simple. 
Still, even for an arbitrary H{q) the precession equations 
have some special properties implied by the 0(3) sym- 
metry that will enable us to find their solutions. In the 
next section we proceed directly to the solution of the 
equations (|17p . In the Appendix we explain the group- 
theoretic foundations of this solution. 



IV. SOLUTION OF THE PRECESSION 
EQUATIONS BY SPINORIAL DECOMPOSITION 

In order to find analytic solutions of the precession 
equations, we express the vector X{q) as a composition of 
a two-component spinor ^(g) and its complex conjugate. 



X{q) = ^^iq)a^{q). 



(18) 



where cr is the vector of Pauli matrices. The precession 
equation (IT71) will be satisfied by the vector X{q) if the 
spinor ^(g) obeys the equation 



E^'iq)=l[cT-n{q)]^iq). 



(19) 



This set of two coupled equations can always be reduced 
to separate second-order equations for the two compo- 
nents of the spinor. This is a crucial simplification — 
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equations of the second-order are much more friendly 
than third-order equations. The technique of lowering 
the order of a differential equation by seeking the so- 
lution in the product form was known before for some 
third-order equations [l^. Here we extend it to a set of 
equations with an underlying symmetry group. 

For our choice of the precession vector we obtain from 
(HH), 



£:V('?) + (l + 9'T»f)V'(9)=0, 



(20) 



where the upper (lower) sign corresponds to the upper 
(lower) component 'ijj{q) of the spinor ^"(9). These equa- 
tions have parabolic cylinder functions as their solutions. 
Of course, the solutions of the second-order equations 
must be chosen in such a way that the spinor ^'(g) is a 
solution of the original equations There are several 
equivalent representations of the solutions of Eqs. (PO)) 
in terms of parabolic cylinder functions (cf., for example, 
[131 ) ■ We have chosen the set of two ortho normal spinors 
in such a way that they contain only two parabolic cylin- 
der functions and their complex conjugates. 



*i = 
where 

F{£,q)=e 



-G{£,q)) 



*2 



G*{£,q) 
F*{£,q)) 



l + i 



G{£,q)^e-^^D_ 



V£ 

l + i 

V£ 



(21) 

(22a) 
(22b) 



The two complex spinors ((2T|) will serve as the building 
blocks in the construction of three real orthonormal vec- 
tors. We shall represent them in a compact way by em- 
ploying the same Dirac matrices (although in a different 
context) as in the formula (j4]). 



1 



where 



$ = 



^-1 

*2 



(23) 



(24) 



and the matrices pk act in the space spanned by spinors 
'I'l and \E'2- For example, Xi = -(\I/|crvI/2 + ^1(t^i)/2. 
Having determined an orthonormal set of vectors that 
solve the precession equations (1171) , we may construct the 
matrix W. It is made of various quadratic expressions 
built of the two functions F and G 



^[F^ + G^] 5R[F2 + G2] -2!3[i^G*] 
-25R[FG] 23[FG] -FF* + GG* 



(25) 



One may check by a direct calculation that this matrix 
is orthonormal W-^W = 1. Owing to the minus sign in 
the definition ([23]) . the determinant of W is equal to 1 
(it describes a pure rotation). The final formula for the 
vector X{q,T) that represents the solution of the initial 
value problem is 

Xiq,T)^Wiq)-[yViq-£r)Y ■X{q-£T,0), (26) 

where the matrices VV are expressed in terms of the 
parabolic cylinder functions F and G by the formula (|25p . 
Note that as compared with the general formula ((T^ we 
could replace the inverse matrix with the transposed ma- 
trix because we used an orthonormal set of solutions. 



V. TIME EVOLUTION OF THE QED VACUUM 



The solution of the initial value problem ([26]) can be 
used to evaluate the overlap Q{q,T) between the DHW 
function evolved from the QED vacuum in time and the 
initial DHW function. This overlap is a simple measure 
of the departure from the initial state. It is given by the 
projection of the vector (|26l) on the vector 



Voiq) = 



1 



q- 



(l,0,g), 



(27) 



that describes the initial vacuum state. In this way we 
obtain the following compact formula for the overlap as 
a scalar product of two vectors. 



where 



Q{q,T) = V{q)-V{q-£T), 



Viq)^W^{q)-Vo{q). 



(28) 



(29) 



The overlap function, denoted by ao in [Ij], was ob- 
tained there by numerical integration of Eqs. ([T]). Our 
analytic formulas reproduce exactly those results but in- 
stead of showing again the plots we shall transform the 
overlap function into a more intuitive characteristic — 
the density of pairs n(q, r) produced by the electric field. 
The relation between these two functions, as shown in [l[ 
by the calculation of the produced current, is very simple 



n(q,r) 



l-Q{q,r) 



(30) 



Obviously, the number of pairs is small when the overlap 
function is close to unity. In Fig. [T] we show the pair 
density as a function of q for different values of r. The 
pair density n(g, r) as a function of q has a plateau in the 
middle region and some transient oscillations near g « 
and q « £t. The length of the plateau region grows with 
T covering more and more of the q axis. The plateau 
region reaches only as far as £t. This is an expression 
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n(q,T) 
0.03 



FIG. 1. (Color online) The pair density n{q, t) as a function of 
the dimensionless momentum q plotted at the values 10, 20, 
30, and 40 of the dimensionless time r. The dimensionless 
field strength £ is set to one. 



of causality because the largest momentum that could 
be acquired in time r is Et. The departures at g = 
and q £t from the plateau values can be explained as 
transient effects. The longer we apply the electric field, 
the less important is their role: they occupy a smaller 
and smaller fraction of the phase space. 

In the limit, when r — oo, the pair density is totally 
dominated by the plateau. The height of the plateau at 
the mid point for q — £t /2 has the value 

(1 - Q(5r/2,r))/2 = (1 - V {Et /2)-V{-£t /2)) 12. 

(31) 

This height can be determined as a function of £, in the 
limit when t — > oo, from the asymptotic expansions of 
the parabolic cylinder functions ([2^ for large values of 
the argument. These expansions are different for q ^ oo 
and q — > — oo due to the Stokes phenomenon. The role of 
the Stokes phenomenon in pair production has been re- 
cently exposed in the context of the WKB approximation 
2l| . The leading terms in the asymptotic expansions are 
11 



-4*2 (g) 



Fi£.-\q\) 



o 



Gi£,-\q\) = e 



where the phases are 

Now, let us note that for |g| ^ 1 the vacuum vector (P7)) 
reduces to (0, 0, sgn(g)). Therefore, the vectors y(±oo) 



(32a) 
(32b) 
(32c) 
(32d) 

(33) 



involve only those elements of the matrix (l25t that form 
its bottom row. This observation leads to a particularly 
simple result for the positive values of q. In this case, 
according to p2ap and (j32bp . we obtain 



y(^) - (0,0,1). 



(34) 



The calculation of V{—oo) is more involved, since the 
asymptotic formulas are more complicated. However, 
because V{oo) has only one nonvanishing component, 
we only need the last component of the vector V(— oo) 
to calculate the scalar product (|28p . This component is 
equal to — |Gp, where we can substitute the asymp- 
totic values given by p2cp and ()32d|) . With the use of 
the following identity for the Euler F function: 



Til + ix)T{l 



IX) 



sinh(7ra;) ' 



we obtain finally 



lim n(fT/2,T) 



oxp(-- 



(35) 



(36) 



It was shown before (cf., Eqs. (69), (72), and (73) in IJ) 
that this value for the number of pairs, after integration 
over all momenta, leads to the Schwinger formula ^15i] . 



VI. CONCLUSIONS 

We have shown once more that the DHW function is 
a powerful tool to analyze quantum electrodynamic pro- 
cesses without recourse to perturbation theory. It is true 
that the evolution equations for this function, in general 
case, are prohibitively complicated. However, the case 
of the time- independent, uniform electric field treated 
here allows for the complete analytic solution of the ini- 
tial value problem. Our derivation of the analytic so- 
lution may also be viewed as a case study to explain a 
new method of solving differential equations with some 
group symmetry. The original equations involved the 
adjoint representation of the 0(3) group. By expressing 
the solutions in terms of the fundamental representation 
(spinors) we were able to simplify the equations. This 
technique will also work for other symmetry groups. 
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Appendix A 

In order to see the underlying group structure, we shall 
rewrite Eq. ([TT]) with the use of the structure constants 
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dj, for the 0(3) group. To show the general structure of 
these equations more clearly, we will change the notation 
as follows: 

^Q'W--c),WWQ''W- (Ai) 

The structure constants for the 0(3) group are given by 
the Levi-Civita symbols e^fc. 

There are two equivalent ways of looking at these equa- 
tions. One may consider the right hand side as the Pois- 
son bracket oiW and G''- Then the equations determine 
a family of canonical transformations. Alternatively, one 
may view them as the evolution equation of a vector X 
in the adjoint representation of the group, i.e. 

-^X{\)=z[S--H{\)]Xi\), (A2) 



where the matrices Si are the spin 1 matrices in the 
Cartesian basis. The second interpretation is closer to 
our procedure. We just replace the evolution equations 
(|A2[) for the vector in the adjoint representation (spin 
1) by the evolution equations for the vector in the fun- 
damental representation (spin 1/2). We only have to 
change the generators S ^ cr/2 but we keep the same 
combination of the generators as determined by H. Af- 
ter solving the equations in the spinorial representation 
we may construct solutions in any other representation 
by taking the appropriate products. 
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